Residual resistance in two-dimensional, microwave driven Hall systems 
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We address the origin of the residual resistance observed in microwave irradiated two-dimensional 
electron gases in a weak magnetic field. We study charge modulations arising from negative photo- 
conductivity and show that dissipative currents are exponentially suppressed. Relating the exponent 
to the temperature dependence of microscopic parameters taken from experiment, we find pseudo 
activated behaviour. In order to obtain these results it is essential to take into account the finite 
range of the Coulomb interaction. 



Recently Mani et al. [l| and Zudov et al. |2( (see also 
Ref. 0) have reported that in ultra-clean Hall systems 
the dissipative resistance is driven to exponentially small 
values by sufficiently strong irradiation with microwaves. 
Theories addressing these unexpected "zero-resistance" 
states are numerous and cover a wide range of ideas 
HUBBJlJi- A calculation by Durst et al. (see 
also Refs. 0,01) revealed that the photo-conductivity os- 
cillates with photon frequency and that it takes negative 
values in regions where, experimentally, the resistance 
is zero. The negative photo-conductivity leads to insta- 
bility towards an inhomogeneous, current carrying state 
0, H| ■ Given such a state the system can sustain a net 
measurement current without voltage drop Q, a result 
which, together with Ref.0, provides an explanation for 
the regions of zero resistance. 

Despite this progress one of the intriguing features 
remains unclear. In the "zero-resistance" regime an 
exponentially small residual resistance remains with a 
temperature dependence, which appears activated like, 
R cx exp(— Tq/T). The measured energies ksTo are sur- 
prisingly large. They exceed the obvious energy scales, 
the cyclotron and microwave energies, by an order of 
magnitude. 

In this letter we particularly address the problem of 
the residual resistance and propose an alternative mech- 
anism to thermal activation. Building mainly on the idea 
of negative local conductivity and the resulting instabil- 
ity, we develop a non-local theory of dissipative conduc- 
tion. We apply methods similar to the ones which have 
been used for studying the Gunn effect [JJj . Within this 
framework we investigate the global spatial structure of 
the electric potential inside the system for given bound- 
ary conditions. We show that dissipative currents are 
suppressed by an exponential factor exp(— A/6), where 
A/2 is the domain size in the charge modulations and b 
is the domain wall width. This is our main result. We 
argue that the residual resistance is suppressed by the 
same factor. Relating b to microscopic parameters mea- 
sured in experiments and taking their observed temper- 
ature dependence into account, we find pseudo activated 
behaviour. In order to obtain these results it is essential 
to use a finite range interaction. 

Our starting point is a phenomenological ansatz for the 



current density, 

j = (a + a p (E)+a H )E-eDySn, (1) 

which contains the electric drift current ooE, the Hall 
current ffijE, and a microwave assisted photo-current 
(Tp(-E)E with E = — the local electric field, and 
E = |E|. D is the diffusion constant, dn is the deviation 
in the electron density from the neutral equilibrium den- 
sity, and e is the electric charge. The unusual transport 
behavior is captured in the non-linear photo-conductivity 
a p {E): in the zero-resistance regime and at vanishing 
electric field <J P (E) is large and negative 0. The exact 
behavior of o~ p (E) depends on many details. However, it 
is clear that, for large electric fields, the photon-current 
saturates and er p (oo) vanishes. Therefore we envisage a 
photo-conductivity that takes a negative minimum value 
at E = and smoothly approaches zero for E — > oo. A 
similar form has been assumed earlier 0. 

Notice that in addition to the currents driven by E 
there is also a diffusion current —eDVSn driven by the 
gradient in Sn. It has to be taken into account, since mi- 
crowave irradiation takes the system out of equilibrium 
and leads to inhomogeneous densities. By writing Eq. 
in this way we have in mind that there is still local equi- 
librium; i.e. conductivities and diffusion coefficient are 
well-defined. Moreover, in the present case they are, to 
a good approximation, independent of Fermi energy and 
temperature 0,0- In this approximation, the relation 
||TJ together with the continuity equation and Poisson's 
equation for the charge density p — eSn provide a closed 
set of non-linear differential equations that determine the 
system states. 

We begin by deriving a single integral equation for the 
electric field in stationary states. To keep things manage- 
able we restrict ourself to one-dimensional (ID) states, 
where the potentials and densities depend only on one 
spatial coordinate, which we take to be x throughout. 
Then, by construction, there is no y-component of E and 
hence no Hall current in the x-direction. In addition, the 
continuity equation requires the x-component of the cur- 
rent to be constant in space, j x — const. — j. Next, we 
integrate Poisson's equation in the form 

4>{x) = [ G{x- x')p{x') dx' , (2) 
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where the reduced ID Green function G(x) is the electric 
potential of a ID unit line-charge p(x, y) = 5(x). We take 
the derivative of this equation with respect to x, integrate 
by parts and use Eq. Q to obtain the non-linear integral 
equation 

E(x) =-^[ G{x - x') [(a + a p (E))E - j] dx' . (3) 

The simplest solutions of this equation are states with 
constant electric field E throughout the system. In this 
case the dissipative current is 

j = [^r + cr + a v {E))E, (4) 
(jo 

where Go = / G(x)dx. This global j — E relation is 
sketched in Fig. QJl. We would like to stress that this 
is not a local current-field relation but strictly holds for 
uniform electric fields only. It follows from a straightfor- 
ward two-dimensional (2D) stability analysis that homo- 
geneous states are stable if and only if \E\ is larger than 
the critical field E defined by D/Go + o- a + a p (E ) = as 
indicated in Fig. ^ The associated Hall current density 
cthEo in y-direction corresponds to the critical current 
density jo in Ref . Q • 

In the limiting case of a short-ranged Green function, 
G(x) — GqS(x), Eq. © becomes local, which allows the 
construction of piecewise constant solutions as depicted 
in Fig.^*. This situation here coincides with the current 
domain picture in Ref. Q. Its foundation on a current- 
field relation like Eq. Q has been pointed out by Begeret 
et al. . Note that only states of vanishing dissipative 
current j and E — ±E can be globally stable (cf. Fig.^l. 
Apart from that, size and configuration of domains of 
critical field ±Eq are largely undetermined. 

The arbitrariness of solutions is an artifact of the zero- 
range interaction. Moreover, in this limit the discussion 
of essential properties of the dissipative currents is not 
possible, as will become clear below. Apart from this, it 
is questionable whether the assumption of a short range 
interaction applies here at all. In fact, the absence of 
metallic gates in the experiments suggests that the 
screening length is large. We therefore set out to con- 
struct global states without having to rely on the as- 
sumption of local interactions. We will also consider 
quasi-stationary states that are weakly time-dependent. 
As above, we restrict ourself to ID states. 

Using the Green function G c of the Coulomb potential 
in Eq. © makes the problem of finding exact solutions 
intractable. However, the model Green function 

G(x) = -W\x\ , (5) 

which is the exact Coulomb propagator of line charges in 
two-dimensional electro-dynamics, allows a surprisingly 
simple construction of non-trivial global solutions [lOj. G 
is a good approximation to G c if one adjusts W such that 




FIG. 1: a) The dissipative current j as function of the electric 
field for homogeneous stationary states according to Eq. @. 
For a given current j only one stable E exists, except for j — 0, 
where one finds two (marginally) stable fields ±Eq. b) Stable 
(solid) and unstable (dashed) potential profiles for a short 
range interaction. The slopes of the stable solution equal 
plus or minus the critical field Eo- The dissipative current j 
of the stable solution vanishes despite the finite potential A(/> 
drop along the system. 

both yield the same force on the relevant length-scale. 
The improvement of choosing G, compared to the delta- 
function approximation, is the finite interaction range. 

Taking the second derivative with respect to x of 
Eq. HJ|, with the original Green function replaced by G, 
reveals that 

-d 2 x (j) = d x E = 2Wp(x) . (6) 

Next, we take the time derivative of this equation, elimi- 
nate dtp using the continuity equation, integrate over x, 
and use Eq. to obtain a differential equation 

dlE - ±d t E = ™ {(cr + a p (E))E - j (t)} (7) 

for E(x,t), the x component of the local, time-dependent 
electric field. The integration constant jo(t) has the 
dimension of a current density. Indeed, for stationary 
states, jo is the constant x-component of the current den- 
sity, as can be verified with help of Eq. Q ■ Actually, in 
the present context an equation very similar to Eq. © 
has recently been discussed by Bergeret et al. 

8]- 

Stationary solutions obey 

d 2 x E - f(E) — J , (8) 

where f(E) = 2W(2a + cr p (E))E/D, and J = 2Wj /D. 
If x is interpreted as time, this is the equation of motion 
of a par ticle with coordinate £ in a force field f(E) — J 

(cf. El). 

For a negative zero-field conductivity — u = ao + °p(0) 
and sufficiently small current J the potential V(E) as- 
sociated with the force in Eq. © has the shape of an 
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FIG. 2: Oscillations of a mechanical particle in the poten- 
tial V(E) translate into periodic modulations of the electric 
potential <j> in the Hall system. For vanishing dissipative cur- 
rent, J = 0, the profile is "flat", while for J > (asymmetric 
oscillations) the potential ascends. 



inverted double- well (Fig. Oscillatory orbits inside 
the well that almost reach the potential maxima at ±Eq 
correspond to periodic field and potential profiles. The 
main differences to the solutions discussed before are the 
rounded edges and well defined domain lengths d_ and 
d + . The rounding width b can be estimated by the fre- 
quency f2 of harmonic oscillations near the minimum of 
the potential, b « 7r/fi. The domain lengths d_ and d + 
are given by the dwelling time of the particle near the left 
or right maximum. They can therefore be estimated by 
d± sa 2 / lu In Eo/q±, where q± is the minimum distance of 
the particle orbit to ±Eq, and ui the frequency of oscil- 
lations in the inverted potential near ±Eq. In parabolic 
approximation of cr p (E) the frequencies u> and £1 differ 
only by a factor \/2, s.t. the estimates for b and d± lead 
to the simple relation 
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Q± 



(9) 



For vanishing current, V(E) is symmetric in E and 
d + = d- . This is the "flat" solution shown in Fig. [5] A 
small but finite dissipative current Sj, however, causes 
an asymmetry in the domain lengths <i_ , d + , given by 
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which in turn gives rise to ascending or descending poten- 
tial profiles (Fig. |2J . The crucial point is that for periods 
A = d- + d+ larger than b an exponentially small dissi- 
pative current j ~ uEq exp(— A/6) comes with a drastic 
change of the potential-profil. 

At first glance, these findings seem to give an explana- 
tion for residual dissipative resistivity: an exponentially 
small dissipative current in x direction with its necessar- 
ily ascending or descending potential profile gives rise to 
a large Hall current (in the y-direction). However, this is 



not quite so simple: the dissipative current of the station- 
ary states flows from the lower to the higher potential, 
as can easily be verified within the mechanical model. 
Therefore the residual resistance ascribed to these states 
is negative. For the same reason these states are dy- 
namically unstable. This can be checked by use of the 
dynamical equation (J7J: A perturbation of the domain 
sizes leads to currents that tend to increase the pertur- 
bation. It is not surprising that the transient currents 
involved in this process are, again, exponentially small, 
oc exp(— X/b). Note that the temporal evolution must 
therefore slow down exponentially with increasing ratio 
X/b. 

Are these transient dissipative currents also responsi- 
ble for the observed residual resistance? To make this 
idea concrete, we analyze in the following one scenario, 
in which we can explicitly relate the transient currents 
to a transient residual resistance of the entire system. 

We consider a symmetric DC bias (Hall voltage) A<j>. 
Knowing that the stationary solutions discussed above 
(cf. Fig. |21) are not physically relevant, it is not clear 
from the start how the potential drop happens across 
the system. We have approached this problem by nu- 
merical simulations of the system's dynamics, governed 
by Eq. , the equation of continuity, and Eq. J2J) using 
the Green function of the Coulomb potential with dif- 
ferent screening lengths s. The lower curves in Fig. |3| 
are the results for zero bias. The resemblance to the 
theoretical profiles with zero current is obvious. Note 
that the profiles keep on developing, albeit exponentially 
slowly with increasing A. The upper curves in Fig. [3] are 
typical potential formations for non-zero bias. By inject- 
ing resp. extracting charges the electric potentials at the 
system boundaries were fixed to given reference voltages. 
Significantly, in the bulk the profile looks the same as 
in the case of vanishing voltage, and the entire potential 
drop happens near the edges over two large edge domains. 
It can also be seen that these domains slowly but steadily 
move inwards. 

It turns out that the transient profile with applied bias 
in Fig. |2| too, can be analyzed in the mechanical analogy, 
when the dynamical term in Eq. (JJJ) is included as an x 
dependent (anti-)friction term. In mechanical terms one 
domain corresponds to a half cycle of oscillation in the 
potential V(E), with the domain length set by the energy 
of the particle. To proceed we use a shock-wave ansatz 
E(x,t) = E(x - vt) HU, with the velocity v piecewise 
constant in x. By means of this the time derivative of 
Eq. becomes a (anti-)friction term —'yd x E(x — vt), 
with coefficient 7 — v/D. This term extracts or supplies 
energy, and allows changes of the domain length at the 
expense of a time dependence in the corresponding po- 
tential profiles. By inspection of Fig. one finds that 
significant friction is needed at the edge domains only. 
In order to switch from a large edge domain length r to 
the smaller bulk period A an exponentially small energy 
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FIG. 3: Numercial results: The electric potential develop- 
ing from an initial random charge density (white noise cor- 
relation) under microwave irradiation at finite (upper) and 
vanishing Hall voltage (lower curve, with a offset of -15 for 
better visibility). The evolution has exponentially slowed 
down (time steps between successive lines increase by a fac- 
tor \/2), however, peaks still tend to merge (structure around 
x « 250). The screening length of s — 65 is larger than the 
period. In response to a finite Hall voltage two large edge 
domains form, which move slowly inwards. Parameters in 
natural units (I = 1, u c = 1, e 2 /l = 1): D — 1.5, St — 0.1, 
u = 0.5, So = 1.0 

of order exp— A/6 has to be removed. One arrives at 
this result with the same ideas that led to Eq. JjJ}. Con- 
sequently, the edge domain must drift with a velocity 
v ~ exp(— A/6) into the system. The additional charges 
needed to move the edge domain imply an exponentially 
small transient current / cx rexp(— X/b). The domain on 
the right edge can be described similarly. Note that for 
large A</> there is Ohmic behavior in the sense that the 
current is proportional to the Hall voltage, because r is 
proportional to A(j>. 

It is apriori not impossible that the observed resisidual 
resistance is in fact transient, as in our example. How- 
ever, this has not been seen yet in experiment. If the 
dissipative current turns out to be persistent, the ID 
treatment presented here cannot provide an answer, and 
a generalization to two dimensions might be necessary. 
Nevertheless, based on the preceding discussion we be- 
lieve that the exponent X/b generally governs the sup- 
pression of the dissipative current, as A and b are the 
only two fundamental lengthscales of the charge modula- 
tion. We expect this to be independent of the particular 
choice of the Green function and screening, provided the 
interaction is long ranged. This is confirmed by numerics 
and further analytical considerations |l3l |. 

We now relate b to experimental parameters and dis- 
cuss below what limits A. For the Coulomb interaction 
we find, by a simple analysis using Eq. J2J, b ~ D/u, 
which is a microscopic length scale of the problem. The 
diffusion coefficient D can be determined by Einstein's 
relation from the dark conductivity. Determining u, the 



dissipative conductivity in the dips, is more difficult, be- 
cause u does not correspond to directly measured quan- 
tities. However, an estimate can be obtained, if one as- 
sumes that the microwave enhancement and reduction 
of the conductivity is symmetric, i.e that it is lowered 
in the dips by the same amount that it is enhanced in 
the peaks (cf. Q). Using the data underlying Fig. 3 in 
Ref. 0, we find b ~ D/u of the order of the magnetic 
length or larger, except for the smallest resistances of the 
l.lfcG minima (see below). Furthermore, the resistivity 
enhancement (solid circles in Fig. 3 in Ref. |2|) is linear in 
the inverse temperature, which suggests the same for the 
reduced resistivities, hence u ~ Cq + C\/T . Recall that 
by our hypothesis the residual resistance is proportional 
to the factor exp (—X/b) ~ exp(— Xu/D). Assuming that 
A is a weak function of T, the C\ /T dependence of u gives 
rise to activated-like behaviour. 

It is clear that the magnetic length is a lower bound to 
b. When this limit is reached at sufficiently low temper- 
atues we expect a drastic change in the T dependence. 
Indeed, this limit seems to be reached in the experiments; 
the estimated value for D/u at the l.lfcG minimum for 
the lowest temperatures in Fig. 3a in Q is about one 
order of magnigtude smaller than the magnetic length. 
A corresponding tendency of leveling off is visible in the 
data. 

Comparing the "activation" temperatures To to 
XCi/D shows that even moderate bulk periods A > b 
can explain the large observed values for Tq. In addition, 
for A and D constant, the linear dependence of C\ on the 
magnetic field strength B, seen in the data Q, implies 
T oc B. 

It remains an open question what determines the do- 
main width A. A finite screening length s is an obvious 
candidate. By incorporating screening into the mechan- 
ical model we found evidence that A can assume certain 
fractions of 4s which, however, can be arbitrarily small 
[l3j |. It is also possible that Landau level depletion, when 
a domain becomes too large, significantly changes the lo- 
cal transport properties to the degree that relation 
simply breaks down. 

A direct measurement of A would be a crucial test for 
the domain picture. Furthermore, experimentally inves- 
tigating correlations between A/6 and the residual resis- 
tance allows a test of our ideas. 

We thank R. G. Mani and M. A. Zudov for supplying 
us with experimental information and data, and K. von 
Klitzing and M. R. Zirnbauer for stimulating discussions. 
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